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Abstract 

Higgs bosons with enhanced coupling to bottom quarks are copiously produced at hadron col- 
liders via bb — > h, where the initial b quarks reside in the proton sea. We revisit the calculation of 
the next-to-leading-order cross section for this process and argue that the appropriate factorization 
scale for the b distribution functions is approximately m/j/4, rather than m^, as had been previously 
assumed. This greatly improves the convergence of the perturbation series, and yields a result with 
mild factorization-scale dependence. We also show that the leading-order calculation of gg — > bbh, 

Mh. 

integrated over the momenta of the final-state particles, is very sensitive to the factorization and 

<D ' 

■ renormalization scales. For scales of order the gg — > bbh cross section is comparable to that 



- 1—1 

X 



of bb — > h, in contrast to the order-of-magnitude discrepancy between these two calculations for 



the scale m^. The result we obtain improves the prospects for Higgs-boson discovery at hadron 
colliders for large values of tan (5. 
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FIG. 1: Leading-order diagram for the production of the Higgs boson via bottom-quark fusion. 



I. INTRODUCTION 

In the standard model, the Higgs boson couples to fermions with strength rrif/v, where 
v = (a/2Gf) -1 / 2 ~ 246 GeV is the vacuum-expectation value of the Higgs field. The Yukawa 
coupling of the Higgs boson to bottom quarks (nib ~ 5 GeV) is thus very weak, leading to 
very small cross sections for associated production of the Higgs boson and bottom quarks 
at the Fermilab Tevatron (y/S = 1.96 TeV pp) and the CERN Large Hadron Collider 
(LHC, = 14 TeV pp) |2|. However, this Yukawa coupling could be considerably enhanced 
in extensions of the standard model with more than one Higgs doublet, thereby increasing 
this production cross section Q]. For example, in a two-Higgs-doublet model, the Yukawa 
coupling of some or all of the Higgs bosons (h°, H°, A , H^) to the bottom quark could be 
enhanced for large values of tan/3 = i^/ui, where v\ is the vacuum expectation value of the 
Higgs doublet that couples to the bottom quark. 

The dominant subprocess for the production of a Higgs boson in association with bottom 
quarks is bottom-quark fusion, bb —>■ h (Fig. Q), 1 where the b quarks reside in the proton 
sea P, [1 0]. This is the leading-order (LO) subprocess for the inclusive production of the 
Higgs boson in association with bottom quarks. Since the bottom-quark sea is generated by 
gluons splitting into nearly-collinear bb pairs, the final state contains two spectator bottom 
quarks that tend to be at low transverse momentum (pr)- 

In contrast, if one requires one bottom quark at high p? from the production process, the 
leading-order subprocess is bg — > bh ja, ui, uj] . This process is particularly promising due to 
the ability to tag the b quark in the final state. The cross section for the production of the 
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gg,qq-^bbh UU\m, IllL ll2L Il3l Il4l Il5l| . 2 Although this process has been the most studied, 



1 We use h to denote a generic Higgs boson. In a two-Higgs-doublet model, h may denote any of the neutral 
Higgs bosons (h°, H°, A°). 

2 qq — > bbh is negligible in comparison with gg — » bbh. 
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FIG. 2: Representative diagrams for associated production of the Higgs boson and two high-py 
bottom quarks: (a) gg — > bbh (8 diagrams); (b) qq — > 66/i (2 diagrams). 

it is likely that bg — *> bh is the more promising, due to its larger cross section. The inclusive 
cross section, bb — > /i, which we study in this paper, is useful when the Higgs boson can be 
identified above backgrounds without the need to detect the accompanying bottom quarks 
that reside in the final state. This subprocess may be useful to discover a Higgs boson for 



large tan/3 in the decay mode h — > t + t at the Tevatron and LHC |l6l. Il7| and h — > fi + fi 
at the LHC 3,Q,Q. It has the advantage of having the largest cross section, since it is 
inclusive of the other two processes. 

Higgs-boson production via bottom-quark fusion was calculated at next-to-leading order 
(NLO) in Refs. There are two puzzling aspects of the results of that calculation: 

• Although the NLO correction is modest, it consists of two independent corrections, of 
order 1/ \n(m h /m b ) and otg, which are both large (and of opposite sign). This suggests 
that the perturbation series in each expansion parameter individually may not be well 
behaved. 

• The cross section at the Tevatron (both LO and NLO) is an order of magnitude larger 
than the cross section obtained by calculating gg — > bbh (Fig. |2J) and integrating over 
the momenta of the final-state particles While gg — ► bbh is not a reliable calcula- 
tion of the total inclusive cross section, since the expansion parameter is as ln(m/j/m&) 
rather than as |3j, the large discrepancy between the two calculations is surprising. 

In this paper we solve both of these puzzles. Implicit in both puzzles is the choice of 
the factorization scale, which had been chosen to be the Higgs-boson mass in Refs. 0, |^. 
Although the choice of the factorization scale in a fixed-order calculation is often regarded 
as arbitrary, we argue that there is a prescription based on physical considerations. Refining 



the discussion of Ref. 



21] , we show that the relevant factorization scale is a fraction of the 



Higgs-boson mass, approximately m h /A. We find that this choice of scale solves both puzzles 



listed above. We thereby present a reliable NLO calculation of the inclusive cross section 
for the production of the Higgs boson in association with bottom quarks. 

In Section Owe determine the relevant factorization scale for Higgs-boson production in 
association with bottom quarks. In Section ITTT1 we present the results of this scale choice. We 
discuss these results in Section II VI and show that they solve the two puzzles listed above. 
Section |V] summarizes the conclusions of our study and suggests further work. 



II. FACTORIZATION SCALE 



The b distribution function, like any parton distribution function, sums (to all orders) 
collinear logarithms that appear at higher orders. Thus, to determine the relevant factor- 
ization scale, we investigate the collinear logarithm that arises at next-to-leading order. 

There are two independent NLO corrections to bb — > h. The first is from initial jrluons, 
bg — > bh (Fig. EJ, which is a correction of order l/ln(mft/m&), as explained in Ref. The 
second is from virtual and real gluon emission (Fig. 0]), which is a correction of order as- 
Since as(m h ) is proportional to 1/ \n(m h / Aq C d), one can regard both of these corrections 
as being of the form of an inverse logarithm. 

The calculations involved are identical to those of Refs. 

a a. 

We keep the effect of 

the bottom-quark mass exactly in our calculations, without approximations. In order to do 
this, we may set the bottom-quark mass to zero in all diagrams in which the bottom quark 
appears as an initial-state parton. 3 This is called the simplified Aivazis-Collins-Olness-Tung 
(ACOT) scheme 22|, |23|, |24j. The only subprocess in which we must keep the bottom- 




s' inroipL 



FIG. 3: Diagrams for the next-to-leading-order correction to bb — > h from initial gluons. This 
correction is of order l/ln(mft/m&). 



3 One may maintain a nonzero bottom-quark mass in these diagrams, but it does not increase the accuracy 
of the calculation. 
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FIG. 4: Diagrams for the next-to-leading-order correction to bb — > h from real and virtual gluon 
emission. This correction is of order as- 

quark mass is gg — > bbh (Fig. EJ), which is a next-to- next-to-leading-order correction of order 
1 / In 2 (mh/rrib). In practice, for rrih ^> m&, it is an excellent approximation to neglect the 
bottom-quark mass in this subprocess; we keep the mass nonzero, nevertheless. 

Let us first investigate the 1/ ]n(mh/mb) correction, from initial gluons. The first diagram 
in Fig. El has a collinear divergence due to the t-channel quark propagator. The hadronic 
differential cross section therefore has the behavior da/dt ~ 1/t in the collinear region. 
The integral over t produces the collinear logarithm. Thus, the upper limit of the collinear 
integration is set by the virtuality, at which the differential cross section begins to 

deviate substantially from the collinear behavior. 

We show in Fig. 03 the hadronic differential cross section times the squared virtuality, 
—t da/dt, versus the virtuality (scaled to the Higgs-boson mass), yf—t/rrih, for a variety of 
Higgs-boson masses at the Tevatron and the LHC. In order to compare the cross sections 
for different Higgs-boson masses at a given collider, we normalize the curves to unity at 
small virtualities. The curves are nearly identical, demonstrating that the differential cross 
section scales with the Higgs-boson mass. At small virtualities the curve is flat, indicating 
the collinear behavior da/dt ~ 1/t. At larger virtualities the cross section is damped. For 
fixed t, the differential cross section is given by []| 

2 1 4,2 



da 
~dt 



1 rs f-±w/s) 
- ds dT]\g(x 1 ,ii F )b{x2,HF)+{xi 

b Jm 2 h -t J\\n{s/S) 



24 I V2 , 



1 m\ + u 1 
s 2 st 
(1) 



which explicitly shows the 1/t behavior for small t. For larger values of —t, the lower limit 
on the s integration, — t, increases and damps the cross section, since the integrand falls 
steeply with increasing s. 

Figure El shows that the virtuality at which the behavior of the differential cross section 
deviates substantially from the collinear behavior is much less than m^. Since the curves vary 
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bg -> bh at 1.96TeVpp 




0.01 0.1 1 2 



bg -> bh at 14 TeV pp 




0.01 0.1 1 2 



FIG. 5: Hadronic differential cross section times the squared virtuality for the subprocess bg — > 6/i 
us. the virtuality (scaled to the Higgs-boson mass) at both the Tevatron (upper plot) and the LHC 
(lower plot). Curves are shown for a variety of Higgs-boson masses, scaled such that they overlap 
at small virtuality. 
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smoothly, there is some ambiguity in defining the virtuality at which the collinear behavior 
ceases. For the sake of discussion, let us define \f—i < rrih/4: as the collinear region. The 
collinear region extends to slightly higher virtualities at the LHC than at the Tevatron, but 
this is small compared with the inherent ambiguity in defining the collinear region. 

The collinear logarithm that is generated at NLO is therefore approximately ln(m/i/4jUp), 
rather than ln(mh/ fip). Thus the factorization scale for the l/ln(m^/m;,) correction should 
be chosen to be of order y,p sa rrih/A in order to sum the collinear logarithm. We will 
examine the consequences of this scale choice in the next section. In order to account for 
the ambiguity in defining the collinear region, we will vary the factorization scale between 
ix f = rrih/8 and [Xp = 

This derivation of the factorization scale is similar to the argument presented in Ref . j^j] , 
where the behavior of the cross section as a function of the transverse momentum of the 
final-state b quark, p T b , is studied. We prefer to instead use the variable which has the 



interpretation of the virtuality of the i-channel bottom-quark propagator. Since we adopt 
the simplified ACOT formalism 22, 0, 3- we are t° set the b mass to zero, which 
makes the discussion of the finite b mass in Ref. [21] moot. Despite these differences, our 



2l| and yields similar results. 



approach is very similar to that of Ref. 

The results for the as correction from virtual and real gluon emission (Fig. HJ) are simi- 
lar. We show in Fig. El the hadronic differential cross section for real gluon emission times 
the squared virtuality, —tda/dt, versus the virtuality (scaled to the Higgs-boson mass), 
a/— t/rrih, for a variety of Higgs-boson masses at the Tevatron and the LHC. A cut of 
— u > (10 GeV) 2 is imposed to regulate the infrared singularity associated with soft-gluon 
emission. The collinear region extends up to a slightly higher virtualities than in the bg — > bh 
subprocess, but the difference is small compared with the inherent ambiguity in defining the 
collinear region. The collinear region again extends to slightly higher virtualities at the LHC 
than at the Tevatron. 



III. RESULTS 

Following Refs. |2l|4|, we fix the renormalization scale of the Yukawa coupling to /xr = rrih- 
It was shown in that paper that the renormalization-scale dependence of the cross section is 
modest, and is reduced at NLO in comparison with LO. Hence we focus on the factorization- 
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bb ->• gh at 1.96TeVpp 




0.01 0.1 1 2 



bb -> p/i at 14 TeV 




FIG. 6: Same as Fig. [SJ but for the subprocess 66 — > gh. A cut of — u > (10 GeV) 2 is imposed to 
regulate the infrared singularity associated with soft-gluon emission. 
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scale dependence of the cross section. 

We show in Fig. [7| the factorization-scale dependence of the inclusive cross section for 
Higgs-boson production via bottom-quark fusion for = 100 GeV at the Tevatron and the 
LHC. The four curves on each plot are described below. The results for heavier Higgs bosons 
are qualitatively similar. While a standard-model Higgs boson of 100 GeV is excluded, the 
lower bounds on the masses of the neutral Higgs bosons h°, A of the minimal supersymmetric 
standard model are about 91 GeV (^j). Even lighter Higgs bosons are possible in a general 
two-Higgs-doublet model j^J. We set tan/3 = 1 throughout. 

Let us first focus on the results at the Tevatron, the upper plot in Fig. The factoriza- 
tion scale covers a wide range, including the canonical choice /zp = used in Refs. 0, 0] 
and fip ~ advocated in the previous section. The curve labeled "LO" is the LO cross 

n 

section calculated with LO parton distribution functions (CTEQ6L1 [27]), which has signifi- 
cant factorization-scale dependence. The curve labeled "LO+1/ In" is the partial NLO cross 
section, including only the l/h^m^/m;,) correction, calculated with NLO parton distribu- 
tion functions (CTEQ6M). At fip = Trih this correction is large and negative, approximately 
—70%. However, at /iF ~ this correction is small, indicating that this is indeed the 

relevant factorization scale for this process. 

The curve labeled "NLO" in Fig. [7| is the full NLO cross section, including both the 
1 / ln(mh / rrib) correction and the as correction, calculated with NLO parton distribution 
functions. At fip = n^h, the as correction is large and positive, nearly canceling the large 
negative 1/ ln(mh/mb) correction. However, at fiF m rrih/4:, the as correction is modest, 
yielding a modest NLO correction. This again indicates that this is the relevant factorization 
scale for this process. The factorization-scale dependence of the NLO cross section is reduced 
in comparison with that of the LO cross section. 

The curve labeled "LO+1/ In +1/ In 2 " in Fig. [TJis the partial NLO cross section [LO plus 
1/ \n(mh / rrib) correction] plus the 1/ \n 2 (mh/mb) correction from the diagrams in Fig. |2j 4 
which is part of the next-to-next-to-leading-order (NNLO) correction. 5 This NNLO correc- 
tion is smallest at fiF m rrih/4:, again indicating that this is the relevant factorization scale 
for this process. This is significant because while it is always possible to find a factoriza- 



4 The correction from qq — > bbh is negligible in comparison with that from gg — ► bbh. 

5 We calculated this curve with NLO parton distribution functions since NNLO parton distribution functions 
are not yet available. 
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FIG. 7: Cross section for Higgs-boson production via bottom-quark fusion vs. the factorization 
scale for = 100 GeV at the Tevatron (upper plot) and the LHC (lower plot). Curves are shown 
for leading-order, next-to-leading order, leading-order plus corrections of order l/ln(m^/mfe), and 
leading-order plus corrections of order 1/ In (mft/mj) and 1/ ln 2 (m/j/mb). 
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tion scale such that the NLO correction is small, it is not guaranteed that this same scale 
will yield a small NNLO correction, unless there is a good motivation for this scale. We 
anticipate that a full NNLO calculation will further support our argument that the relevant 
factorization scale for this process is fip ~ Wjj/4. 

The results at the LHC, the lower plot in Fig. are qualitatively similar to those at 
the Tevatron. We argued in the previous section that the relevant factorization scale for 
this process at the LHC is slightly higher than at the Tevatron. The partial NLO cross 
section [LO plus 1/ ' \n{mh/m^) correction] crosses the LO cross section at a slightly higher 
factorization scale at the LHC than at the Tevatron, consistent with this argument. 

We present in Tables ITHITTlthe NLO cross sections for Higgs-boson production via bottom- 
quark fusion at the Tevatron (both y/S =1.8 and 1.96 TeV) and the LHC, using fip — ^/i/4. 
These cross sections differ from those of Refs. ^, 0] in part due to the improved choice of 
factorization scale, and also in part due to a bug in the CTEQ4M j3] computer code that 
affected the gluon and b distribution functions used in that paper. Since there is some 
ambiguity in defining the collinear region, we vary the factorization scale between twice 
and one-half its central value, and consider this an uncertainty in our calculation. This 
corresponds to the first uncertainty listed in Tables U MlIIl We also vary the renormalization 
scale between and 2m^, and report this as the second uncertainty in Tables ITV lIIII This 

uncertainty is considerably less than that associated with the factorization scale. 

There are two additional sources of uncertainty in our calculation. The uncertainty in 
the 6-quark MS mass, m^m?,) = 4.24 ± 0.11 |29j, yields an uncertainty in the Yukawa 
coupling (evaluated at [ir = m^). This corresponds to the third uncertainty in Tables U tillH 
which is the same for all Higgs-boson masses and machine energies. The fourth uncertainty 
corresponds to the uncertainty in t he p arton distribution functions, which we evaluated 
using the method described in Refs. {2 7^ [3^] . The four sources of uncertainty are combined 
in quadrature and reported as an absolute uncertainty on the NLO cross section. The 
same exercise is performed for the LO cross section, where we report only the combined 
uncertainty. 
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TABLE I: Leading-order (LO) and next-to-leading-order (NLO) cross sections (pb) for Higgs-boson 
production via bottom-quark fusion at the Tevatron (\/S = 1.8 TeV pp). The LO cross sections are 
computed using CTEQ6L1 parton distribution functions |27j ] and 1-loop evolution of the bottom- 
quark Yukawa coupling, ybip-n)- The NLO cross sections are computed using CTEQ6M parton 
distribution functions and 2- loop evolution of the Yukawa coupling and as (fin)- The four sources 
of uncertainty in the NLO cross sections are also listed. The factorization scale is fip = m h/4, 
and is varied between fj,p = m^/8 (upper uncertainty) and fj,p = (lower uncertainty). The 

renormalization scale is hr = m^, and is varied between [ir = (upper uncertainty) and 

Hp = 2rrih (lower uncertainty). The uncertainty in the Yukawa coupling stems from the uncertainty 
in the b mass, m& = 4.24 ± 0.11 GeV. The final uncertainty is due to the uncertainty in the parton 
distribution functions. These four uncertainties are combined in quadrature and reported as an 
absolute uncertainty in the NLO cross section. The combined uncertainty in the LO cross section 
is also given. 
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TABLE II: Same as Table ID but for y/S = 1.96 TeV. 



m h (GeV) a LO (pb) a NLO (pb) 5a NLO (5fi F ,5fi R ,5y b , PDF) (%) 
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-3 


8.87 


+1.84 
-1.65 


xlO" 


-3 


-12 
+13 


±3.8 


±6.4 


+14 
-12 


125 


5.72 


+2.42 
-2.63 


xl(T 


-3 


7.34 


+1.55 
-1.36 


xlO - 


-3 


-11 
+13 


±3.7 


±6.4 


+15 
-13 


130 


4.76 


+1.96 
-2.14 


xl(T 


-3 


6.11 


+1.32 
-1.14 


xlO" 


-3 


-11 
+12 


±3.7 


±6.4 


+16 
-14 


140 


3.34 


+1.32 
-1.43 


xlO" 


-3 


4.28 


+0.97 
-0.81 


xlO~ 


-3 


-9.3 
+12 


±3.5 


±6.4 


+18 
-15 


150 


2.37 


+0.91 
-0.98 


xl(T 


-3 


3.04 


+0.73 
-0.59 


xl(T 


-3 


-8.3 
+11 


±3.4 


±6.4 


+20 

-16 


160 


1.70 


+0.65 
-0.69 


xlO - 


-3 


2.19 


+0.56 
-0.45 


xl(T 


-3 


-7.5 
+10 


±3.3 


±6.4 


+22 
-18 


180 


9.10 


+3.47 
-3.52 


xlO" 


-4 


1.18 


+0.34 
-0.26 


xlO" 


-3 


-6.2 
+9.6 


±3.1 


±6.4 


+26 
-20 


200 


5.05 


+2.00 
-1.92 


xlO" 


-4 


6.66 


+2.19 
-1.64 


xlO" 


-4 


-5.2 
+8.9 


±2.9 


±6.4 


+31 

-23 


250 


1.30 


+0.62 
-0.51 


xl(T 


-4 


1.81 


+0.81 
-0.57 


xlO - 


-4 


-3.6 
+7.6 


±2.7 


±6.4 


+43 
-30 


300 


3.83 


+2.30 
-1.64 


xl(T 


-5 


5.69 


+3.31 
-2.17 


XlO" 


-5 


-2.9 
+6.7 


±2.5 


±6.4 


+57 

-37 


400 


4.21 


+3.88 
-2.26 


xlO" 


-6 


7.48 


+6.85 
-3.87 


XlO" 


-6 


-2.6 
+5.2 


±2.2 


±6.4 


+91 
-51 


500 


5.62 


+7.59 
-3.69 


xlO - 


-7 


1.23 


+1.66 
-0.80 


xlO" 


-6 


-2.4 
+4.6 


±2.1 


±6.4 


+134 
-64 



IV. DISCUSSION 

The first puzzle listed in Section H] is solved by choosing the factorization scale appropri- 
ately. As we showed in the previous section, for /xp w the 1/ \a{rnh/m%) correction is 
small and the as correction is modest. Thus perturbation theory in each expansion param- 
eter individually is well behaved. Furthermore, we evaluated the 1/ In 2 (m^/mb) correction 
as well, and found that it is vanishingly small at the relevant factorization scale, providing 
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TABLE III: Same as Table E but for the LHC (y/S = 14 TeV pp). 



m h (GeV) (Tlo (pb) ctnlo (pb) 5ct N lo(<W, <W, PDF) (%) 



60 


2. 


.87 


+3.44 
-2.32 


Xl0° 




4.69 


+1.70 
-2.35 


Xl0° 




-49 
+34 


±6.8 


±6.4 


+5.9 
-7.3 


70 


2. 


.04 


+2.05 
-1.51 


Xl0° 




3.23 


+0.99 
-1.35 


XlO° 




-40 
+29 


±6.3 


±6. 


.4 


+5.5 
-7.0 


80 


1. 


.48 


+1.29 
-1.01 


Xl0° 




2.29 


+0.61 
-0.82 


XlO° 




-34 
+25 


±5.8 


±6. 


.4 


+5.2 
-6.7 


90 


1. 


.09 


+0.85 
-0.70 


Xl0° 




1.65 


+0.39 
-0.52 


xlO° 




-30 
+22 


±5.6 


±6. 


.4 


+5.0 
-6.4 


100 


8. 


.20 


+5.79 
-4.91 


Xl0~ 


l 


1.22 


+0.26 
-0.34 


xlO° 




-26 
+19 


±5.4 


±6. 


.4 


+4.7 
-6.2 


105 


7. 


.15 


+4.83 
-4.16 


Xl0~ 


l 


1.06 


+0.22 
-0.28 


xlO° 




-25 
+18 


±5.3 


±6. 


.4 


+4.6 
-6.1 


110 


6. 


.26 


+4.06 
-3.54 


Xl0~ 


l 


9.19 


+1.83 
-2.34 


xlO~ 


i 


-23 
+18 


±5.2 


±6. 


.4 


+4.6 
-6.0 


115 


5. 


.50 


+3.44 
-3.03 


xl0~ 


l 


8.04 


+1.55 
-1.96 


xlO~ 


i 


-22 
+17 


±5.1 


±6. 


.4 


+4.5 
-5.9 


120 


4.85 


+2.92 
-2.61 


xl0~ 


l 


7.05 


+1.31 
-1.64 


xlO~ 


i 


-21 
+16 


±5.0 


±6. 


.4 


+4.4 
-5.8 


125 


4. 


.29 


+2.50 
-2.25 


xl0~ 


l 


6.21 


+1.12 
-1.39 


xlO~ 


i 


-20 
+16 


±4.9 


±6. 


.4 


+4.3 
-5.7 


130 


3. 


.81 


+2.15 
-1.96 


xl0~ 


l 


5.48 


+0.96 
-1.18 


xlO~ 


i 


-19 
+15 


±4.9 


±6. 


.4 


+4.3 
-5.6 


140 


3. 


.03 


+1.61 
-1.49 


xl0~ 


l 


4.32 


+0.71 
-0.87 


xlO~ 


i 


-18 
+14 


±4.7 


±6. 


.4 


+4.1 

-5.5 


150 


2. 


.44 


+1.22 
-1.15 


xlO~ 


l 


3.45 


+0.54 
-0.65 


xlO" 


i 


-16 
+13 


±4.6 


±6. 


.4 


+4.0 
-5.3 


160 


1. 


.98 


+0.94 
-0.90 


xlO~ 


l 


2.78 


+0.42 
-0.49 


xlO~ 


i 


-15 
+12 


±4.5 


±6. 


.4 


+3.9 
-5.2 


180 


1. 


,34 


+0.58 
-0.57 


xlO~ 


l 


1.86 


+0.26 
-0.30 


xlO~ 


i 


-13 
+11 


±4.2 


±6. 


.4 


+3.8 
-4.9 


200 


9. 


.31 


+3.77 
-3.73 


xlO~ 


2 


1.29 


+0.17 
-0.19 


xlO~ 


i 


-12 
+10 


±4.1 


±6. 


.4 


+3.7 
-4.7 


250 


4. 


.19 


+1.45 
-1.48 


xlO~ 


2 


5.69 


+0.66 
-0.70 


xlO~ 


2 


-8.9 
+8.3 


±3.8 


±6. 


.4 


+3.7 
-4.4 


300 


2. 


.11 


+0.64 
-0.67 


xlO~ 


2 


2.83 


+0.31 
-0.31 


xlO~ 


■2 


-7.1 
+7.2 


±3.6 


±6. 


.4 


+3.9 
-4.3 


400 


6. 


.66 


+1.69 
-1.81 


xlO~ 


3 


8.81 


+0.90 
-0.88 


xlO~ 


3 


-4.9 
+5.6 


±3.3 


±6. 


.4 


+4.9 
-5.2 


500 


2. 


.56 


+0.58 
-0.62 


xlO~ 


3 


3.37 


+0.36 
-0.34 


xlO~ 


3 


-3.6 
+5.2 


±3.1 


±6. 


.4 


+6.5 
-6.5 


600 


1. 


.13 


+0.24 
-0.25 


xlO~ 


3 


1.47 


+0.17 
-0.16 


xlO~ 


3 


-2.8 
+4.5 


±2.9 


±6. 


.4 


+8.6 
-8.0 


700 


5. 


.42 


+1.11 
-1.16 


xlO~ 


4 


7.12 


+0.95 
-0.85 


xlO~ 


4 


-2.3 
+3.8 


±2.8 


±6. 


.4 


+11 

-10 


800 


2. 


.80 


+0.59 
-0.59 


xlO~ 


4 


3.69 


+0.57 
-0.49 


xlO~ 


4 


-1.9 
+3.5 


±2.7 


±6. 


.4 


+13 
-11 


900 


1. 


.53 


+0.34 
-0.32 


xlO~ 


4 


2.03 


+0.36 
-0.30 


xlO~ 


4 


-1.5 
+3.1 


±2.6 


±6. 


.4 


+16 
-13 


1000 


8. 


.73 


+2.05 
-1.86 


xlO - 


5 


1.17 


+0.23 
-0.19 


xlO~ 


4 


-1.2 
+2.8 


±2.5 


±6. 


.4 


+19 
-15 
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further evidence that the perturbative series in l/ln(mft/m&) is wel 
terminates at this order, while the series in as extends to all orders 



aehaved. This series 



The second puzzle is also solved by a consideration of the choice of scales, both factor- 

Finn 

ization and renormalization. In Refs. |3|, y,|20j], all scales were chosen to be the Higgs-boson 
mass, /ip = /ir = m/j. The NLO cross section at the Tevatron is nearly a factor of ten 
greater than that obtained by calculating gg — > bbh and integrating over the momenta of 
the final-state particles. However, this factor is much less at lower scales, mostly because the 
cross section for gg — > bbh is very scale dependent, and increases significantly at lower scales. 
In contrast, the NLO cross section for bb — > h has mild scale dependence, and decreases by 
only about 25% for \i F rs m h /4:. 

We have established that the relevant factorization scale for bb —>■ h is /ip ~ rrih/A. It 
is likely that the relevant factorization scale for gg — > bbh is also much less than rrih, as 
well as the renormalization scale of as- We have made no attempt to establish the relevant 
renormalization scale for the Yukawa coupling in bb — > h, and we have found that our NLO 
calculation is insensitive to this scale. However, the size of the as correction is less for 
smaller renormalization scales, which suggests that the relevant renormalization scale may 
be less than rrih- Let us adopt this ansatz, although we do not have a rigorous justification 
for it, in contrast to our derivation of the factorization scale. 

As a specific example, we evaluate the cross section for gg — > bbh with fip = fiR = rrih 
and \if = = rrih/4: for rrih = 100 GeV at the Tevatron, and compare with our NLO 
calculation of bb — > h. The results are listed in Table IIVI The order of magnitude difference 
between our NLO calculation and gg — > bbh when the scale is rrih is reduced to about a 
factor of two for scales close to rrih/ 4. A factor of two can easily be accounted for by the 
fact that our calculation sums collinear logarithms to all orders, while gg — > bbh produces 
only the LO collinear logarithm. This solves the second puzzle listed in the introduction. 

It is the desire to sum collinear logarithms that leads one to use bb — > h as the LO 
subprocess for the inclusive production of the Higgs boson. The corrections are of order 
1/ hi(mh/mb) and as, and we have seen that they are modest for the appropriate choice 
of factorization scale. If one were to use gg — > bbh as the LO subprocess, the expansion 
parameter would be a5ln(m/,/m(,), and perturbation theory would be poorly behaved. Our 
calculation gives the most accurate and reliable cross section for the inclusive production of 
the Higgs boson because it sums these collinear logarithms to all orders. 
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TABLE IV: Cross sections (fb) for bb -> h at NLO and gg -> at LO for m h = 100 GeV at the 
Tevatron, for two choices of the common factorization and renormalization scales. The final column 
gives the ratio of the cross sections. The ratio is nearly an order of magnitude for /zp = [ir = m>h, 
but is only about a factor of two for fip = [ir = 













Scales 


a{bb -► h) 


a{gg bbh) 


a(bb 


h)/a(gg -> 6Wi) 


fJ-F = fJ-R = m h 


26.6 fb 


3.1 fb 




8.5 


fJ-F = Vr = mh/4: 


20.8 fb 


9.2 fb 




2.3 



It is suggested in Refs. that the calculation of bb — ► fa may overestimate the 

inclusive cross section, due to crude approximations inherent in the kinematics, which give 
rise to larg e bottom-quark mass and phase-space effects. However, the ACOT formalism 
22i l23l l24j makes no approximations in either the kinematics or the b mass; we maintained 
the effect of the b mass exactly. We find no evidence for any inconsistency in the ACOT 
formalism. Rather, we find that the LO calculation of gg — > bbh, integrated over the 
momenta of the final-state particles, underestimates the inclusive cross section when the 
factorization and renormalization scales are chosen to be m^. 



V. CONCLUSIONS 



We have revisited the next-to-leading-order (NLO) calculation of Higgs-boson production 

nn 

via bottom-quark fusion and solved the two puzzles associated with that calculation (3, |4| . 
We showed that the appropriate factorization scale for this process is fip ~ m /i/4, rather 
than m/j, as had been previously assumed. This greatly improves the convergence of the 
perturbation series, which was mediocre for fip = m^. The resulting cross section has mild 
factorization-scale dependence, and small renormalizat ion-scale dependence. It is the most 
reliable calculation of the inclusive cross section for Higgs-boson production in association 
with b quarks. 

To support our arguments, we calculated one of the next-to-next-to-leading (NNLO) 
corrections (associated with the diagrams in Fig. |2J), and showed that it is vanishingly small 
for jip ~ TUh/4:. The ingredients exist to calculate the full NNLO cross section, using the 
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results from Refs. [?], Q|- This should yield a cross section with small factorization- and 
renormalization-scale dependence. It will also provide an additional check of our choice of 
factorization scale. 

The other puzzle we solved also involves the choice of scales, both factorization and 
renormalization. The inclusive cross section for Higgs-boson production in association with 
bottom quarks may be approximated by gg — > bbh, integrated over the momenta of the final- 
state particles. This yields a result an order of magnitude less than the NLO calculation of 
bb — > h. However, this result is very scale dependent, since it is based on a leading-order 
calculation. Choosing scales of order rather than m^, we find that the cross section 

is comparable to that of bb — > h (see Table HV|) . The NLO calculation of gg — > bbh might 
support these observations. However, that calculation is not as accurate as one based on 
bb — > h for rrih ^> nib, since the latter sums collinear logarithms to all orders in perturbation 
theory. 

Let us review the existing calculations of Higgs-boson production in association with 
bottom quarks. The relevant calculation depends upon the final state that is desired. For 
the inclusive cross section, the relevant leading-order (LO) subprocess is bb — > h (Fig.HJ. The 



NLO cross section was calculated in Refs. |3|, and updated in this paper. The cross section 
for the production of the Higgs boson accompanied by one high-transverse-momentum (jpx) 
bottom quark is obtained at LO from the subprocess bg — ► bh (Fig. EJ), which is calculated 
at NLO in Ref. |2j. This process is particularly promising due to the ability to tag the b 
quark in the final state. Finally, the cross section for the production of the Higgs boson 
accompanied by two high-pr b quarks is obtained at LO from the subprocesses gg,qq — > bbh 
(Fig. |2J). This process has been calculated only at LO thus far, but the ingredients exist to 



provide the NLO cross section |35, 



37, 



38, 



39|. Although this process has been the most 



studied, it is likely that bg — > bh is the more promising, due to its larger cross section. The 
inclusive cross section, bb — > h, obtained in this paper, is useful when the Higgs boson can 
be identified above backgrounds without the need to tag b quarks, such as h — > t + t~, ji + . 
It has the advantage of having the largest cross section, since it is inclusive of the other two 
processes. 
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